Application of the lent particle method 
to Poisson driven SDE's 
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Abstract 

We apply the Dirichlet forms version of Malliavin calculus to stochastic differen- 
tial equations with jumps. As in the continuous case this weakens significantly the 
assumptions on the coefficients of the SDE. In spite of the use of the Dirichlet forms 
theory, this approach brings also an important simplification which was not available 
nor visible previously : an explicit formula giving the carre du champ matrix, i.e. the 
Malliavin matrix. Following this formula a new procedure appears, called the lent 
particle method which shortens the computations both theoretically and in concrete 
examples. In this paper which uses the construction done in [7] we restrict ourselves 
to the existence of densities, smoothness will be studied separately. 

AMS 2000 subject classifications: Primary 60G57, 60H05 ; secondary 60J45,60G51 
Keywords: stochastic differential equation, Poisson functional, Dirichlet Form, energy 
image density, Levy processes, gradient, carre du champ 

1 Introduction 

During the last twenty years a significant development of the theory of Dirichlet forms 
occured in the direction of improving regularity results or lightening hypotheses of Malli- 
avin calculus ^3]. With respect to the Malliavin analysis on Wiener space, what brings 
the Dirichlet forms approach is threefold: a) The arguments hold under only Lipschitz 
hypotheses, e.g. for regularity of solutions of stochastic differential equations cf [6], this 
is due to the celebrated property that contractions operate on Dirichlet forms and the 
Emile Picard iterated scheme may be performed under the Dirichlet norm, b) A general 
criterion exists, the energy image density property (EID), proved on the Wiener space for 
the Ornstein-Uhlenbeck form, and in several other cases (but still a conjecture in general 
since 1986 cf [5]), which provides an efficient tool for obtaining existence of densities in 
stochastic calculus, c) Dirichlet forms are easy to construct in the infinite dimensional 
frameworks encountered in probability theory (cf [6] Chap.V) and this yields a theory of 
errors propagation, especially for finance and physics cf [3], but also for numerical analysis 
of PDE's and SPDE's cf [3l]. 

Extensions of Malliavin calculus to the case of stochastic differential equations with 
jumps have been soon proposed and gave rise to an extensive literature. The approach is 
either dealing with local operators acting on the size of the jumps (cf [2] [TT] [24] etc.) or 
acting on the instants of the jumps (cf [lO] [l2]) or based on the Fock space representation 
of the Poisson space and finite difference operators (cf [26] [27] [18] etc.). In all cases the 
arguments are somewhat intricate. 
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We have obtained results which simphfy highly the approach with local operators cf 
[7]. Based on Dirichlet forms on the general Poisson space they gather the advantages of 
Dirichlet forms methods and simplicity of use. It may be summarized in the following way: 
in order to calculate the Malliavin matrix, we add a particle to the system, compute the 
gradient of the functional on this particle, and take back the particle before integrating by 
the Poisson measure. The main formula is 



T[F] = [ e~{-f[e+F])dN. 
Jx 



where 7 is the carre du champ operator on the state space, T the carre du champ operator 
on the Poisson space and e"*", e~ are the operations of adding and cancelling a particle. 

Let us present this method, called the lent particle method, on a simple example. Let Y 
be a real Levy process with absolutely continuous Levy measure v = kdx, thus such that 
1 + AYs / almost surely. We equip the space M \ {0} with the Dirichlet form on L'^{y) 
with carre du champ operator 

7[n](.x) =x2'u'^(x)l^|^|<i}. 

The link of Y with the associated random Poisson measure N with intensity dt x a 
is V/i G L^^^(M"'"), Jq h{s) dYs = f l[Q^t]{s)h{s)xN{dsdx) where N is the compensated 
measure N — dt x a. 

Let us consider the Doleans-Dade exponential 

£{Y)t = e''^Yl{l + AYs)e'^^^. (1) 

s ^ t 

In order to study the regularity of the pair (Yf, £{Y)t), we proceed as follows: 

1. We add a particle (a,y) i.e. a jump to Y at time a with size y what gives 

{£iY)t) = e^*+^ + An)e^^= (1 + y)e-y = £{YUl + y). 

s ^ t 

2. We compute ^[e+£(Y)t] = {£{Y)tfyH^^y^^iy 

3. We take back the particle before integrating w.r.t. N : 

e-^[e+£{Y),] = {£iYUl + y)-')\\^^^iy 

and 



my)t] = E i£iY)t{i + AY^)-TAY^^^y^^^,y 



a < t 

As easily seen, by a similar computation, the matrix T of the pair {Yt,£{Yt)) is given by 



r_V-f 1 £{YUl + AY^)-^ \ 



a ^ t 
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Hence the density of the pair {Yt,S{Yt)) under hypotheses yielding (EID) lies on the con- 
dition 



where JT denotes the jump times of Y with size less that ^ between and t. This is 

guaranteed if Y has an infinite Levy measure. 

Hence, summarizing our hypotheses (see Lemma [1] below) 

Let Y be a real Levy process whose Levy measure is infinite, absolutely continuous w.r.t. 
Lebesgue measure and whose density dominates a positive continuous function near 0, then 
the pair {Yf,£{Y)t) possesses a density on M^. 

The aim of the present article is to apply the lent particle method to stochastic differ- 
ential equations with jumps. 

We shall begin by explaining the method (Sections 2.1 and 2.2). The only difference 
with our treatment in [7] is that we work here on a product probability space in order to 
be able to put also a Brownian motion or other semi-martingales in the studied SDE. With 
respect to the note [4| we have introduced a clearer new notation, as in [7] the operator e~ 
is shared from the integration by N . 

The SDE we consider is described in Section 2.3 with the assumptions done. 

The argument dealing to the Malliavin matrix is concentrated in Section 2.6 thanks to 
the lent particle method. 

We give some examples in Section 3 of SDE driven by Levy processes when the Levy 
measure is underdimensioned and the diffusion matrix is degenerated. This improves 
known results on the Levy area (Section 3.2). We end by lightening assumptions on exis- 
tence of density for a McKean-Vlasov type non-linear SDE (Section 3.3) and for stable-like 
processes (Section 3.4). 

2 Notations and hypotheses. 

2.1 Dirichlet structures and Poisson measures. 

Let [X, X, u, d, 7) be a local symmetric Dirichlet structure which admits a carre du champ 
operator i.e. {X,X,v) is a measured space, v is a-finite and the bilinear form 



is a local Dirichlet form with domain d C L'^{v) and carre du champ operator 7 (see 
Bouleau-Hirsch [6j, Chap. I). We assume that for all x G X, {x} belongs to X and that v 
is diffuse {v{{x}) = Vx). The structure (X, A', z^, d, 7) is called the bottom structure. 

We are given N a Poisson random measure on [0, -|-oo[xX with intensity dt x v[du) 
defined on the probability space (rii,^i,Pi) where Oi is the configuration space, Ai the 
a-field generated by and Pi the law of A^. We set iV = iV - x z^. 

We suppose that the bottom structure and A'^ satisfy the following hypothesis : 

Hypothesis (HO). The structure (X, A", d, 7) with generator (a, P(a)) is such that 
there exists a subspace H of ^(a) f] -L^(z^), dense in L^{i>) n L'^{u) and such that V/ e 





3 



Hypothesis (HO) implies what we call the bottom core hypothesis in [7] and denote (BC). 
It is a technical condition due to the fact that the carre du champ takes its values in 
and we need a set of test functions for which it has its values in L^. This condition is not 
so restrictive. In the case of a Poisson measure induced by a real Levy process whose Levy 
measure is absolutely continuous w.r.t. the Lebesgue measure, the following Lemma gives 
a way to fulfill it: 

Lemma 1. Let r G N*, {X,X) = {W , B{W)) and v = kdx where k is non-negative and 
Borelian. We are given ^ = i^r an W^^ -valued and symmetric Borel function. 

We assume that there exist an open set O <Z W and a function ip continuous on O and 
null onW\0 such that 

1. u{dO) = 0, 

2. k > on O v-a.e. and is locally bounded on O 

3. ^ is locally bounded and locally elliptic on O in the sense that for any compact subset 
K in O, there exists positive constants ck and Ck such that 

T 

Vx = (xi, • • • G i^, ^ \ii,3{x)\ ^ Ck and ^ iij{x)xiXj ^ cx|x|^, 
I. k ip > u-a.e. on O 

5. for all i,j G {1, • • • ,r}, belongs to Hl^^{0). 

We denote by H the subspace of functions f G L^{i^) H L^{v) such that the restriction of 
f to O belongs to C^{0) (i.e. C°° with compact support in O). Then, the bilinear form 
defined by 




is closable in L'^{y). Its closure, (d,e), is a local Dirichlet form on L'^{v) which admits a 
carre du champ 7. Moreover, it satisfies hypothesis (HO) and property (EID) i.e. for any 
d and for any -valued function U whose components are in the domain of the form 

C/4(det7[C/,C/*]) -z^] < A'^ 

where det denotes the determinant and X'^ the Lebesgue measure on iW^ ,B{W^)) . 

Proof. First of all, since u{dO) = 0, it is clear that H is dense in L'^{u) n L^{v). We have 

V/, g G H, e{f, g) = - [ a{f){x)g{x) v{dx), 

Jw 

where a(/) = lo{x){k{x))~^ J2i j 9j{£,i,ji/jdif). Since k ^ tp on O, it is clear that for any 
compact set K C O, k~^ is bounded on K. Prom this, it is clear that a{f) belongs to 
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LP'{v). So, a is an operator with dense domain, symmetric and —a is non-negative hence it 
is closable. As a consequence, the bihnear form (if, e) is also closable and it is clear that 
its closure is a Dirichlet form on Lp'{v) and that the carre du champ operator 7 is given by 



r 



ip{x) 



V/ G H, 7(/)(x) = kji^)dJix)d,fix) 



k{x) 



with the convention g = 0. 

This ensures that hypothesis (HO) is satisfied. 

The (EID) property can be proved using standard arguments (see Theorem 2 in □ 

Hypothesis (HO) is satisfied in several other cases not so simple to describe (cf [15] 
Section 3.1). 

When the aim is to obtain density for M'^- valued random variables, d > 1, and only in 
that case, we have to suppose additional conditions : 
Hypothesis (HI): The structure {X,X,h',d,'y) satisfies (EID). 

Hypothesis (H2): X admits a partition of the form: X = -BlJdJfcS^fc) where for all k, 
Ak £ X with v^Ak) < +00 and v{B) = 0, in such a way that for any k £ N* may be 
defined a local Dirichlet structure with carre du champ: 



with V/ e d, f\A^ e dfc and 7[/]|Afe = 7/[/[aJ- 

Hypothesis (H3): Any finite product of structures Sk satisfies (EID). 

The need for hypotheses (HI) to (H3) lies in the argument followed in [7] to prove 
(EID) on the Poisson space. As mentionned above, no case is known where (EID) fails for 
a local Dirichlet structure with carre du champ. In the classical applications hypotheses 
(HO) to (H3) are fulfilled. 

We consider also another probability space {^,2,^.2,^2) on which an M"- valued semi- 
martingale Z = {Z^, - ■ ■ ,Z^) is defined, re G N*. We adopt the following assumption on 
the bracket of Z and on the total variation of its finite variation part. It is satisfied if both 
are dominated by the Lebesgue measure, uniformly w.r.t. w: 
Assumption on Z: 

There exists a positive constant C such that for any square integrable ]R"-valued pre- 
dictable process h: 



The presence of this semimartingale and of the following framework is due to the form 
of the SDE that we will study below cf equation ([8)). This form will be discussed in §2.3. 
We shall work on the product probability space: 




{n, A, P) = {Qi X !^2, ^1 8) ^2, Fi X F2). 
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As in [7], starting from the Dirichlet structure on the bottom space {X, X, v) we construct a 
Dirichlet form on L?'{Q.i,Ai,¥i) and then by considering the product of this Dirichlet struc- 
ture with the trivial one on L^(02, IP2), we obtain a Dirichlet structure (fi, P, B, £, V) 
with domain B C L'^{VL,A,¥) and carre du champ operator T . As we assume (HO) to (H3) 
we know that it satisfies (EID) (cf [7] or Song [32]). 

Finally we denote by {At) ^ the natural filtration of the Poisson random measure N on 
[0,+oo[xX 

2.2 Expression for the gradient and the carre du champ operator 

Following [7] §3.2.2, we are given an auxiliary probability space {R, TZ, p) such that the 
dimension of the vector space L'^{R,Tl, p) is infinite and we construct a random Poisson 
measure NQ p on [0, +oo[xX x R with compensator dtx i/x p such that if = - £{ai,ui) 
then N Q p = £{ai,ui,ri) where (r^) is a sequence of i.i.d. random variables independent 
of A^ whose common law is p and defined on some probability space {Cl, A, P) so that NQp 
is defined on the product probability space: {Q,A,F) x {{l,A,F). 

We assume that the Hilbert space d is separable so that the bottom Dirichlet structure 
admits a gradient operator, -D, and we choose a version of it with values in the space 
Ll{R,n,p)={gGL^{R,n,pyjj^gir)p{dr)=0}. 
We denote it by b. 

Let us recall some important properties: 

• Vti G d, Du = v}' e L^{X X R,X (g,n,u X p). 

• Vti G d, Jj^ Wv}' {■ , r) p{dr) = 'y[u]. 

• (chain rule in dim 1) if F : M — > M is Lipschitz then Vn G d, {F o = {F' o u)u^. 

• (chain rule in dim d) if F is (continuously differentiable) and Lipschitz from M"' 
into M then 

d 

Vn = (m, • • • , n,) G d^ (F o u)' = J^iFl ° 

1=1 

Finally, although not necessary, we assume for simplicity that constants belong to dioc (see 
Bouleau-Hirsch [6] Chap. I Definition 7.1.3.) and that 

1 G dioc which implies 7[1] = and l'' = 0. (3) 

We now introduce the creation and annihilation operators and £~: 

V(t, u) G [0, +oo[xX, Vu?i G ill, = 

{{t,u)GSUpp 101 } • 

In a natural way, we extend these operators on by setting ii w = {wi,W2)- 

4,")^'^) = (4,u)(^i)'^2) and e^^„)(t/^) = {e^^^^-^{wi),W2), 
and then to the functionals by 

e'^ H{w,t,u) = H{£f .w,t,u) and e~ H{w,t,u) = H{e7. .w,t,u). 



6 



We now recall the main Theorem of [7] which gives an explicit formula for a gradient of 
the upper structure {Q,A,F,0,£,r). 

Let us introduce some notations. The space Vq was introduced in hypothesis (HO). We 
denote B the completion of Pq ® -^^^([0, +00 [, dt) d with respect to the norm 



POO p \ 2 r 

\H\\d = (E/ / e"{-f[H]){w,t,u)N{dt,du)) + E / / {e~\H\){w,t,u)'n{t,u)N{dt,du) 
Jo Jx J Jo Jx 



00 



E/ / j[H]{w,t,u)i^{du)dt] +E / / \H\{w,t,u)r]{t,u)u{du)dt, 
Jo Jx J Jo Jx 

where 77 is a fixed positive function in L^(M+ x X,dt x dv). 

Finally we denote by Pat the measure Pat = W{dw)Nw{dt,du). One has to remember that 
the image of P x x by e"*" is nothing but Pat whose image hy e" is 'F x u x dt (see 
Lemma 13 in [7]). 

Theorem 2. The Dirichlet form (B,i5) admits a gradient operator that we denote by jj and 
given by the following formula: 

/•+00 r- 

VFgD, F« = / / e"((e+F)^)(iiV0p G l2(P X P). (4) 

Jo JxxR 

Formula ^ is justified by the following decomposition: 

FeO e+F-FGl ''M^ {{e+ F)') e lKf^x p) F^ G L^FxI 

where each operator is continuous on the range of the preceding one and where Lq(Fn x p) 
is the closed set of elements G in L'^{Fn x p) such that J^Gdp = PAr-a.e. 
Moreover, we have for all F gB) 



r+oo r 

T[F]=E{F^f= / / e- F])dN, (5) 
Jo Jx 



f +00 

/O JX 

where E denotes the expectation with respect to probability P. 

Proof. This is a slight modification of Theorem 17 in [7]. □ 
Let us recall without proof some properties of this structure which are quite natural. 

Proposition 3. If h € L'^{R+,dt) ® d, then N{h) = J^°° h{t,u)N{ds,du) belongs to 
D and 

T[N{h)] = / / ^[h{t,-)]{u)N{dt,du). (6) 

JO JX 

(N{h)) = / h\t,u,r)NQp{dt,du,dr). (7) 
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Remark 1. Theorem 1 gives a method for obtaining r[F] for F G B or F G D", then with 
the hypotheses giving (EID) it suffices to prove det T[F] > P-a.s. to assert that F has 
a density on M". Let us mention a stronger condition which may be also usefull in some 
apphcations. By the following lemma that we leave to the reader 

Lemma 4. Let he random symmetric positive matrices and fi{da) a random positive 
measure. Then {det f MaH{da) = 0} C {/ detMafJ-ida) = 0}, 

it is enough to have / det e~ {'j[e'^ F])dN > F-a.s. hence enough that det e^{'y[e~^F]) 
be > P/vT-a.e. We obtain finally, by lemma 13 of [7], that a sufficient condition for the 
density of F is det -yle^F] > P x z/ x dt-a.e. (or equivalently that the components of the 
vector (e^F)^ be P x z> x dt-a.e. linearly independent in ). 



2.3 The SDE we consider. 

Let d eW, we consider the following SDE : 



Xt = x+ [ I c{s,Xg-,u)N{ds,du) + [ a{s,Xs-)dZs (8) 
Jo Jx Jo 



X 



below denoted (R). 



d ^ ■ Tn>+ TOrf X X ^ R'^ and cj : M+ X ^ M'^^" satisfy the set of hypotheses 



Comment. With respect to the classical form of the SDE's related to Markov processes 
with jumps, as e.g. Ikeda-Watanabe [T7| Chap IV §9, we put the Brownian part in the 
semi-martingale Z. Let us emphasize that the Malliavin calculus that we construct, does 
not concern Z but only N. In fact Z could be replaced by a more general random measure 
with hypotheses assuring existence and uniqueness of the solution. In most applications 
we have in mind there is no Brownian motion at all, since otherwise this induces strong 
regularity properties and classical Malliavin calculus applies. 

Assumption (R): For simplicity, we fix all along this article a finite terminal time T > 0. 
We assume that N and Z are as explained in §2.1. We suppose equation ([8]) is such that : 

1. There exists rj G L?'{X,v) such that: 

a) for all t G [0, T] and u G X, c(t, is differentiable with continuous derivative and 

Vn G X, sup \Dxc{t,x,u)\ ^ fy(ti), 

tG[0,T],xGK'* 

b) \l{t,u) G [0,r] X f/, |c(t,0,u)| ^ 

c) for all t G [0, T] and x G M'', c(t, x, •) G d and 

sup -i[c{t,x,-)]{u) ^ 

te[o,T],xeR'* 

d) for all t G [0, T], all x G M"^ and u G X, the matrix I + -D^c(t, X, u) is invertible and 



sup 

t£[0,T],x& 



{I + DxC{t,X,u)) ^ i^7]{u). 



8 



2. For all t G [0,T] , cr(t, •) is differentiable with continuous derivative and 

sup \Dx(j(t,x)\ < +00. 

t€[0,T],xGR'i 

3. As a consequence of hypotheses 1. and 2. above, it is well known that equation ([H) 
admits a unique solution X such that E[sup^g[o,T] l^iP] < +00. We suppose that for all 
t G [0,T], the matrix {I + '^^^^ Dxa.j{t, Xf.-)AZ-l) is invertible and its inverse is bounded 
by a deterministic constant uniformly with respect to t G [0,T]. 

Remark 2. Assumption (R.3) is satisfied if for example one assumes that there exists a 
constant a > such that for all t G [0, t], all x G X and all j G {1, • • • , n} 

\AZj\ ^ a and l-DxC. j{t,x)\ ^ — . 

na 

But it may be verified also without these inequalities. For example when Z is a Levy 
process, because of the independence of the jumps to the past, the invertibility of (/ + 
Yl^=i DxO'-j{t, Xf-)AZll) is guaranteed in the case d = nhy the assumption that the Levy 
measure has a density on and it remains only to verify the bound of the inverse. 

2.4 Spaces of processes 

We denote by V the predictable sigma-field on [0, T] x Q and we define the following sets 
of processes: 

• TC : the set of real valued processes {Xt)te[o,T]: defined on {Q,A,F), which belong to 
L2([0,T] X f]). 

• Ti-p : the set of predictable processes in TC. 



: the set of real valued processes (-H^t)te[o,T]i which belong to L ([0,T];]D)) i.e. 
such that 



\H\\^^^ = E[ / \Ht\^dt] + / £{Ht)dt < +CX). 
Jo Jo 



\Ht\^dt\ + / 
'0 Jo 

• "^B.-p : the subvector space of predictable processes in 

• Ti-ii®d,v ■ the set of real valued processes H defined on [0, T] x Q x X which are 
predictable and belong to L^([0,T];B d) i.e. such that 

\Mnn^^v=^ir I \Ht?^{du)dt\+ f [ £{Ht{.,u))u{du)dt+E[r e{Ht)dt]<+cx,. 
Jo Jx Jo Jx Jo 

We define H.^ p to be the set of elementary processes in Tin^-p of the form 

m— 1 

Gt{w)=Y,Fi{^)Mu,U+,]{t), 



i=0 
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where m G N*, ^ to ^ " " " *m =^ T and for all Fj S B and is -measurable. 

We also consider "Hp^^j the set of elementary processes in 7{o(^d,v' H belongs to Ti-jj,^^ -p 

if and only if 

m— 1 

Ht{w,u) = Fi{w)l]t^^t,+i]it)9iiu), 

where for alH G {0, • • • , m — 1}, Fj G O and is ^(.-measurable and gi G d. 
The proof of the following lemma is straightforward: 

Lemma 5. ^ is dense in TLn^v "'^'^ '^nxgjd "P dense in 'Ho(gid,v- 
Let Y £ Ti-D^-p, c and a satisfying the set of assumptions (R), then 

{t,w,u) — > c{t,Yt{w),u) G 7^D(g)d,•p 

{t,w) — >a{t,Yt{w)) £nB,P- 

Notation : We shall consider R'^-valued processes, so W^^n^,..., will denote the spaces of 
M*^- valued processes such that each coordinate belongs respectively to 7i,7in,- ■ ■, equipped 
with the standard norm of the product topology. 

The above spaces yield the following results allowing to perform the Emile Picard iteration 
procedure with respect to the Dirichlet norm : 

Proposition 6. Let H G Wu^d.-p and G G TC^ -p, then: 

1. The process 

VtG[0,T], Xt = [ [ H{s,w,u)N{ds,du) 
Jo Jx 

is a square integrable martingale which belongs to Tin and such that the process X~ = 
(X(-)(g[o,T] belongs to 7io,v- The gradient operator satisfies for all t G [0,r].- 

Xl{w,w)= [ [ H\s,w,u,w)dN{ds,du)+ [ [ H\s,w,u,r)NQp{ds,du,dr). 

Jo Jx Jo JxxR 

(9) 

Moreover 

VtG[0,r], IIX^IIb ^V2\\H\\n„^^^^ (10) 
II X W-Hj, ^ V2f II H \\HB(i,d,v (11) 

and 

II \\Hb,v =^ II ^ IIWra®d.P • (12) 

2. The process 

VtG [0,r], Yt= [ G{s,w)dZ, 
Jo 

is a square integrable semimartingale which belongs to Tio, Y^ = (^-)tG[o,T] belongs 
to TCo^-p and 

rt 

yte[0,T], Yliw,w) = G^{s,w,w)dZs. (13) 

^0 
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We also have the following estimates: 

Vi G [0,r], II Yt lb ^C\\G ll^n^ and \\ Y,- ||b ^ C || G H^^.^, 

II Y \\n, ^CVt\\G ||«n^ and \\ ||r,,,^ ^ CVT \\ G ||^n ^ . 
Proof. 1. Assume first that H G W{^g,jj-p, 



m— 1 



H{w,s,u) = Fi{w)l]t^^t,+^^{s)gi{u). 

1=0 

Then, X is a square integrable martingale and we have: 

yt G [0,T],Xt = YFiN{l]t^^t,u+iM]-9i)- 

i 

So, Xt belongs to B and thanks to the functional calculus : 



(n, r)N p{ds, du, dr) 



H^{w, s, u, w)dN{ds, du) + 



]tiAt,ti+iAt] JXxR 

H^{w, s, u, r)N p{ds, du, dr). 



JX 



JXxR 



This yields: 



nxt) 



E 



Jx 



+ E 



JXxR 



H\s,u,'w)\'^ v{du)ds 
I £{H{s,u)u{du))ds +E[ [ e{H{s,-))ds] 

\J0 JX JO 



\H {s,u,r)\^u{du)p{dr)ds 



so we obtain inequalities ([TO]) and (fTTl) in this case and then for general H G 'Hv,(^a,v by 
density. 

As both X and X'^ are stochastic integrals w.r.t. random measure, we deduce that 
\\ms^t,s<tXt exists in -L^(0,P), \\m.s^t,s<tXl exists in L^(rixJ7,PxP) so that, by standard 
arguments, \\m.s^t,s<t exists in D and as we have 

Vs<t, II X, lie ^V2\\H\\n„^^^^, 

what yields this part of the proposition. 

2. For the second part, consider first G G (Wjj-p)"" 

m—l 

Gt{w) = Y,F^{w)l^U,U^^]{t). 
i=0 
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We have 

m—l 



4 = 

m— 1 
1=0 

Thanks to the bounds ([2]), we easily obtain 

E[\Ytf] ^ CE [ \Gsfds, 
Jo 

and 

1 fi2 ft ft " 

£{Yt) = -EE[|y/|2] ^ —EE \Gi\'ds = J2£{Gl)ds. 

It is now easy to conclude, using a density argument, similarly to the proof of the first 
part. □ 



2.5 Preliminary results on the solution of equation 

Proposition 7. The equation ([8]) admits a unique solution X in 7i^. Moreover, the 
gradient of X satisfies: 



t r rt 



Xl = / / D^c{s,Xs^,u) ■ xl^N{ds,du) + / c\s,Xs-.,u,r)N Q p{ds,du,dr) 

Jo Ju Jo JxxR 



+ [' D^a{s,Xs-)-Xl_dZs 
Jo 

Proof. We define inductively a sequence {X'") of M'^- valued semimartingales by X^ = x 
and 

Vr G N, Vt G [0,r], = x + / / c(s, X^. ,u)N{ds, du) + [ a{s,Xl^)dZs. 

Jo Jx Jo 

As a consequence of Proposition [6l it is clear that for all r, X"^ belongs to Ti^ and that we 
have Vt E [0, T] 

= [ [ D^c{s,Xl_,u)-Xl'lN{ds,du)+ [ [ c\s,Xl_,u,r)N Q p{ds,du,dr) 

Jo Ju Jo JuxR 

+ f DMs,x:_)-x:^dZs. 

Jo 

This is the iteration procedure due to Emile Picard and it is well-known that 

lim ^[ sup \Xt-Xl\'^]=0. (14) 

r^+oo te[0,T] 
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Moreover, thanks to the hypotheses we made on the coefficients, it is easily seen that there 
exists a constant k such that for all r G N* and all t £ [0, T] 



EE 



'■+l,tt|2 
t I 



^ K 1 + 



SO that by induction we deduce 

Vr G N, Vi G [0,r], EE 



IX. 



EE 



r,h2 



ds 



< Ke 



Kt 



Hence, the sequence {X^') is bounded in which is an Hilbert space. Therefore, there is 
a sequence of convex combinations of X"^ which converges to a process Y G 7i^. But, by 
(fT4l) we a priori know that X*" tends to X in L'^{[0,T];R'^) so that Y is nothing but X. 
This proves that X belongs to and the relation satisfied by the gradient is consequence 
of relations ((9|) and □ 

We can now explicit a formula for the carre du champ operator of X^, using the linear 
equation satisfied by X". The obtained formula (Theorem [9] below) is known for long time 
(cf [2]) but established here under much weaker regularity assumptions similar to those of 
[llj . Let us emphasize that we obtain this formula without the intensity of the Poisson 
measure being the Lebesgue measure as supposed by these authors. 

Let us define the M"^^ '^-valued process Us by 



Then the following 



tidxd 



dUs = Y,D,a.,,{s,Xs-)dZi. 
i=i 

valued process is the derivative of the flow generated by X: 



Dxc{s,Xs-,u)Ks_N{ds,du) + / dUsKg- 
lo Jx Jo 

Under our hypotheses, for all i ^ 0, the matrix Kt is invertible as a consequence of the 
following proposition which extends classical formulas about linear equations (e.g. [28] 
Chap V §9 Thm 52). 

Proposition 8. Let Hit be a d x d-matrix semimartingale such that L + AT,t is invertible 
\ft a.s. Let Kt be the solution of 

ft 



Kt = I + 



dTi^K, 



(15) 



a) Then Kt is invertible and its inverse Kt satisfies 

Kt = I- [ V^,_(AS,)2(/ + AS,)-i+ [ A',d<S^S^>, . 

b) The d x 1-solution ofVt = Va + f^^ dTisVs- is given by Vt = KtKaVa- 

c) Let Rt be a d X 1-semimartingale. The solution of St = Rt + J^dTisSg- is given by 

St = Ro + Kt[ f Ks-dRs - V Ks-AHs{I + A^s^^ ARs - f Ksd < >,]. 
Jo Jo 



(16) 
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Proof. By Ito's formula the solutions to ([T5]) and (fTBl) satisfy 
dKtKt = d^tKt-Kt^ - Kt^Kt-dllt 

+Kt^Kt^d < S'^ >t - < dT.^K^.K^dT.'' >t 

this equation is Lipschitz with respect to the unknown KK and is verified by the identity 
matrix, what gives the first assertion of the proposition. The remaining is similar. □ 

Here Kt = {Kt)~^ satisfies: 

Kt = I- I [ KsAl + Dxc{s,Xs--,u)y^D^c{s,Xs-,u)N{ds,du) 
Jo Jx 

-[ Ks-dUs+y^Ks-{AUsf{I + AUs)-^+[ Ksd<U'',U^>s. 

■^0 s<t -^0 



2.6 Obtaining the Malliavin matrix thanks to the lent particle method. 
Theorem 9. For all t £ [0, T], 

T[Xt] = Kt f I KMs, X,_, ■)]!<: N{ds,du)K;, 
Jo Jx 

where for any matrix D, D* denotes its transpose. 

Proof. We apply the lent particle method: 

Let {a,u) £ [0,T] x X. We put X^^"'"^ = eJ^^-jXt. We have, 

X^"'") = x+ r [ cis,xl''_'''\u')Nids,du')+ r a{s,xl''_^''^)dZs + c{a,xl"^''\u) 
Jo Jx ^ Jo ° 

+ / [ c{s,X^^'''\u')N{ds,du')+ [ 

J]a,t] Jx * J]a,t] * 

Let us remark that xj:"'^^ = Xt if t < a so that, taking the gradient with respect to the 
variable u, we obtain: 



J]a,t] Jx 



+ [ DMs,x^;_^^^).{x^;_'''^)'dZs. 

Let us now introduce the process K^'^''^^ = e~l^^^^{Kt) which satisfies the following SDE: 



'0 Jx Jo 
and its inverse K^^'^^ = {K^°^'^^)~^ . Then, using the flow property, (Prop ^h), we have: 
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(Let us note that at this stage Remark 1 of Section 2.2 could be used to get a sufficient 
condition of density of Xt.) 

Keeping in mind that the measures P x and Pat are mutually singular, let us emphasize 
that this result is an equality P x h'{du) x p-a.e. (cf Thm 1 above or [7j Prop 18). Now, 
we calculate the carre du champ and then we take back the particle: 

which is an equality PAr-a.e. (cf Thm 1 or [7] Prop 18). 
Finally integrating with respect to N we get 

yt^O,T[Xt] = Kt f I Ks7[c{s,X,-,-)]{u)K:N{ds,du)K; 
Jo Jx 

which ends the proof. □ 
Remark 3. We could also write T[Xt] as 

Kt f j K,- {I+D^c{s, X,- , n))-S[c(s, X,- , ■)]{u){I+{D.,c{s, X,- , u)Y)-^Kl^ N{ds, du)Kl 
Jo Jx 

Remark 4. With the assumptions (R) that we have taken on the SDE ([8]), Xf is in D 
and T[Xt\ given by the above formulae is in L^(P). Now, let us recall that there exists 
a powerfull Borelian localization procedure in any local Dirichlet structure with carre du 
champ (cf [6] Chap. I §7.1). In practice, assumptions (R) may be lightened in such a 
way that significant positive quantities be only finite almost everywhere, then the above 
formulae are still true with Xt £ O/oc and T[Xt] finite a.s. and property (EID) applies as 
well. 



3 Some applications. 
3.1 The regular case 

This is the case where we assume that X is a topological space and that coefficient c is 
regular with respect to the jumps size. More precisely, we have the following 

Proposition 10. Assume that X is a topological space, that the intensity measure ds x u 
of N is such that v has an infinite mass near some point uq in X . If the matrix {s, y, u) — > 
7[c(s,y, •)](u) is continuous on a neighborhood of {0,x,uq) and invertible at {0,x,uq), then 
the solution Xt of JS]) has a density for all t g]0,T]. 

Proof. Let us fix t g]0, T]. As u has infinite mass near uq, as X is right continuous and 
7[c] continuous, admits almost surely a jump at time s S]0,t] with size u £ X such that 
^[c{s, Xg- , ■)]{u) is invertible. As a consequence, 

T[Xt] ^ KtKsl[c{s,X,-,-)]{u)K:K:, 

for the relation order in the set of non-negative symmetric matrixes. As T[Xt] dominates 
an invertible matrix, it is also invertible and this permits to conclude. □ 

Now, the method can yield existence of density for Xt even when the matrix 7[c(s, x, •)] 
is everywhere singular. A example of such a situation was given in [7] §5.3. 
We now turn out to study other degenerated examples. 
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3.2 Levy's stochastic area. 

This example will show that the method can detect densities even when both the matrix 
7[c] is non invertible and the Levy measure singular. 

Let X(t) = {Xi{t), X2{t)) be a Levy process with values in with Levy measure a. 
We suppose that the hypotheses of the method are fulfilled, we shall explicit this later on. 

Let us consider for the moment a general gradient on the bottom space : 

f = /(6 + /26 

where = and ^i, ^2 are functions defined on x which satisfy: r)p{dr) = 

lRC2{-,r)p{dr) = 0, J^^f{xi,X2,r)p{dr) = q;ii(xi,X2), ffi^i{xi,X2,r)C2{xi,X2,r)p{dr) = 
ai2(xi,X2), ^^il{xi,X2)p{dr) = 022(2:1, 2:2), so that 

l[f] = aufi + 2ai2f[f2 + "22 /2^- 

Let be ^ ^ 

V ={Xi{t),X2{t)J Xi{s^)dX2{s)- [ X2{s^)dX,{s)). 
Jo Jo 

We have for < a < t and x = (xi, X2) G M^, 

^{a,x)^ = ^ + {xi,X2,Xi{a-)x2 + xi{X2{t) - X2{a)) - X2{a.)xi - X2{Xi{t) - Xi{a)) 

= V + (rEi,X2,xi(X2(t) - 2X2(0)) - X2(Xi(t) - 2X1 (a))) 
because e~^V is defined P x vxda-a.e. and x da is diffuse, so 

(e+vf = (ei,6,6(^2(t) - 2X2(0)) -6(^1 (i) - 2X1 (a))) 

and 

(ail ai2 Aaii — Bai2 

ai2 a22 Aai2 - Ba22 

Aaii — Bai2 Aa\2 — Ba22 A'^an — 2ABai2 + -6^022 

denoting A = (X2(t) - 2X2(0)) et B = {Xi{t) - 2Xi{a)). 

The operator e~ gives a functional defined PAr-a.e. so that for example 

^("«,xi,x2)(^(*)) = ^(*) - ¥N{dadxidx2)-&.e. 

This yields 

e~A = X2{t) - AX2(a) - 2X2 (a_) let us denote it A 
e'B = Xi{t) - AXi{a) - 2Xi(a_) let us denote it B 

and eventually 

/ aii(AX^) ai2(AX„) iaii(AX„) - ^ai2(AX„) 

r[^]=I] ~ a22(AX„) iai2(AX„)-Sa22(AX,) 

~ ~ A'^au{AXa)-2ABai2{AXa) + B'^a22{^Xa) 

the symbol ~ denoting the symmetry of the matrix. 
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3.2.1. First case. 



Let us consider the case ai2 = 0. We are in this case if v possesses a density satisfying our 
hypotheses which are fulfilled as soon as we assume those of Lemma [U and under these 
assumptions, with same notation, we may choose ol\\ = 022 = {{xl + ^^Dnfy)- have 




r[v]= y lAx^i^ y \+\Ax, 




a < t 

Hence X has a density if the dimension of the vector space spanned by 

I , I 1 \ ,aGJT 

X2{t) - AX2(a) - 2X2{a-) J \ Xi{t) - AXi(a) - 2Xi(a-) 

is equal to 3, where JT = {a G [0, t], AX^ G O}. 

Let us suppose that y{0) = +00. Then one of the projections of O on the axes has an 
infinite mass. Let us suppose it is that of Xi. 

The process Xi{t) — /S.Xi{a) — 2Xi{a—) = Xi{t) — Xi[a) — Xi{a—) cannot remain 
constant for a G JT since the Levy measure of Xi is infinite, and therefore JT is dense in 
M+ and each point of JT is the limit of an increasing sequence (a^) of points of JT such 
that |AXi(Qfc)| — > 0. The constancy of Xi{t) — AXi(a) — 2Xi{a—) would imply that of 
—2Xi{a—) and there will be no jumps. Thus, V has a density if the Levy measure, v , of 
X satisfies hypotheses of Lemma\^ and ^{0) = +00. 

Example 1. Let us take the Levy measure of (Xi,X2) expressed in polar coordinates as 

K^i/9,d0)=5Wde.l]o,i[(/')- 

P 

with g locally bounded and such that it dominates a continuous and positive function 
near 0. Then V = {Xi{t),X2{t), Xi{s-)dX2{s) - X2{s-)dXi{s)) has a density (and 
condition (0.4) of [9] or of [2^ propl.l are not fulfilled). 



3.2.2. Second case. 

Let us suppose ^2 = A(xi, X2)Ci- We are in this case if the measure u is carried by a graph 
in and image of a measure on M. For instance if X2 is taken to be [Xi] the Levy measure 
is carried by the graph X2 = x\. 

Then a\2 = Aan et 022 = A^an. We have 

/I A i-A^ \ 

r[1/] = ^ aii(AX„) A A2 Ai-A2^ . 

a^t \A-\B \A-\^B {A-\Bf I 
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where A is taken on the jumps : X{AXa). If an ^ 0, V has a density as soon as 

dim £ I I A j , a G JT | = 3 (17) 

\\a-xbJ J 

with A-XB = -X2{a^) + X{AX{a))Xi{a-) + X2{t)-X2{a)-X{AX{a)){Xi{t)-Xi{a)). 
In order to study condition (flTl) let us reason on the set 





f I \ 








1 , Jt\ 




\ A-XB j 





A = {oj : dim C 

There exist a, 6, c (depending on oj) such that Va G JT 

a + bX + c{A- XB) = 0. 

If the Levy measure of X is infinite, JT is dense in M+, each point of JT is hmit of an 
increasing sequence (a^) of points of JT such that |AX(afc)| — > 0. If the function A goes 
to zero at zero, the process a + c{X2{t) — 2X2(0—)) vanishes on JT . If the Levy measure 
of X2 is also infinite, ^2(0— ) cannot remain constant when a varies, hence c = a = h = 

Thus, choosing an = x\f\l (or an = Xl^^^ if one assumes hypotheses of Lemma\^, 
if X tends to zero at zero, if the Levy measure of X2 is infinite, and such that there exists 
a bottom structure (M^\{0}, a, d, 7) allowing (BC) and (EID) on the upper space, then 
V has a density. 

Example 2. This applies to V = (Xi(t), [Xi]t, Xi{s^)d[Xi]{s) - Jl[Xi]{s^)dXi{s)). 

The Levy measure of {Xi, [Xi]) is carried by the curve X2 = x\. We have A(a;i,a;2) = 
2x1. We arrive to the sufficient condition : V has a density as soon as the Levy measure 
of Xi satisfies hypotheses of Lemma\^with ^{O) = +00. 



3.3 McKean-Vlasov type equation driven by a Levy process. 

The following nonlinear stochastic differential equation 

{ Xt = Xo + /o a{Xs- ,Ps)dYs te [0, T] 
\ Vs S [0, T], Pg is the probability law of Xs 

where F is a Levy process with values in M"^, independent of Xq, and a : M'^ x'P(M'^) M'^^'' 
where V(R'^) denotes the set of probability measures on M.^, generalizes the McKean-Vlasov 
model. It has been studied by Jourdain, Meleard and Woyczynski [20] who proved, by 
a fixed point argument, that equation (fTSl) admits a solution as soon as a is Lipschitz 
continuous on M'^ x V{M.^) equipped with the product of the canonical metrics on and 
a modified Wasserstein metrics on probability measures. 

When y is a one-dimensional Levy process and k = 1, these authors obtained the 
existence of a density for Xt using a Malliavin calculus in the Bichteler-Gravereaux-Jacod 
spirit under the assumptions that a does not vanish, admits two bounded derivatives with 
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respect to the first variable, and the Levy measure of Y dominates an absolutely continuous 
measure with C^-density and infinite mass, and additional technical conditions. 

We would like to illustrate the lent particle method by simplifying their proof and 
lightening the hypotheses. 

The clever remark — evident after a while of reflection — used by these authors, and 
usefull for us too, is that as soon existence for equation (fTSl) has been proved, it may be 
considered for Malliavin calculus as an equation of the form 



Xt = Xo+ [ a{Xs-,s) dYs 
Jo 



which is a particular case of our present study. 
Let us proceed with the following hypotheses : 

(i) a is n Lip with respect to the first variable uniformly in s and 

sup\{I + D^a)-^{x,t)\ ^ r] 

t,x 

(ii) the Levy measure of Y is such that a Dirichlet structure may be chosen such that 
(HO) and (EID) be fulfilled on the Poisson space (we shall detail this assumption later on). 

By the lent particle method we obtain 



r[Xt] = Kt 



A'„a(X„_,a)7[i,j1(Ay,)a*(X,_,a)< 



k; (19) 



where JT is the random set of jump times of Y before t, 7 is the carre du champ of the 
bottom structure, j is the identity map on W^, Kt is the solution of 

Kt = I+ I dZ, Ks- 



Jo 

where dZg = Yl'i=i Dxa.i{Xs~, s) dY^, and K the inverse of K. 

As a consequence of formula (fT9]l we can conclude that Xt possesses a density on 
under the following hypotheses : 

1") the Levy measure of Y satisfies hypotheses of Lemma [U with ^{0) = +00. Then 
we may choose the operator 7 to be 

'^^^) i=i i=i 

and the identity map j belongs to d and 7[j,i*](x) = (see [7] for a weaker 

assumption and the proof of (EID) on the Poisson space). 

2°) a satisfies (i), is continuous with respect to the second variable at 0, and such that 
the matrix aa*{XQ,0) is invertible. 
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3.4 Stable-like processes. 

The passage between a Levy kernel x, dy) with which is expressed the generator of a 
Markov process with jumps, with standard notation 



^^=1 . X (20) 

+ / \ f{x + y)- fix) - l|y|<i V dif{x)yi v{t, X, dy) 

iK^{o} V y 

to the Poisson random measure N{dt, du) to be used for the SDE able to yield the Markov 
process as solution 

dXt = a{t,Xt)dBt + b{t,Xt)dt + j c{t, Xt-,u)N{dt,du) 

is theoretically always possible thanks to a result of El Karoui-Lepeltier [13]. But this 
general procedure yields for c{t, x, u) a function with few regularity (as the theorem allowing 
to simulate any probability law on M*^ thanks to a random variable defined on [0, 1] equipped 
with the Lebesgue measure, cf [8] Chap. I Thm A. 3.1). 

The study of the correspondence between u{t,x^dy) and the pair {c{t,x,u), N{dt,du)) 
has been deepened by Tsuchiya [33] in order to find conditions on vit, x, dy) such that 
a function c{t, x, u) may be found verifying the Lipschitz hypotheses guaranteeing the 
existence of a solution to the SDE, hence of the Markov process with generator (f2Ql) . 

He applies this study to the case of so-called stable-like processes of order a{x) intro- 
duced in dimension one by Bass [1], whose generator may be symbolically written 

a(x) 

L = -(-A)— (21) 

and he obtains the existence of the Markov process for a Lipschitz and such that ^ a{x) < 
2. (see also [E] Chap. 4). 

This example is pushed further by Hiraba [16] who studies the density of the corre- 
sponding Markov semi-group using the Malliavin calculus in combining the approaches of 
[2] and [21]. He obtains the existence of a density under the hypothesis that a{x) is 
bounded with bounded derivatives and < Ai ^ a{x) ^ A2 < 2 (see also Negoro [25] for 
an analytic proof of this result under C°° assumptions). 

The method of Dirichlet forms allows several improvements of this subject. First about 
the correspondence between vit, x, dy) and the pair (c(t, x, u), N{dt, du)) it is not necessary 
to deal with a Poisson measure whose intensity be the Lebesgue measure but only that it 
carries a Dirichlet form satisfying hypotheses (HO) to (H3). The choice of the operator 7 
is also flexible. 

Second the existence of densities may be performed under weaker assumptions on the 
function c{t, x, u). 

Let us write explicitly this example. The operator 

Af{x) = [ {fix + y)- fix) - D,fix).yl\y\^,)Ciaix))—^ 
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represents the symbol ([2T]) if Ae''^'^ = e''^'^{—\u\"^^^) and this is realized if the function C 
is such that 



what gives (cf [T6] ) 

C p = sm — — ■ 

As these authors let us except the large jumps for the sake of simplicity. We have to obtain 
the operator 

A'' fix) = [ {fix + y)- fix) - D,J.y)Kix, dy) 

J\y\<uo 

with 

dr 



!• poo 

Kix,A) = Ciaix)) da 1a\{o}M 



isd- 

thanks to an SDE of the form 

ft 



Xt = x + / ciXs-,a, z)Nids,da,dz). 

Jo JSd-ixR+ 

Here 5^^"^ is the unit sphere in and da the area measure. If we choose the intensity of A'' 
to be dtdadz on 1R+ x S"^'^ xM_|_ and the function c(x, a, z) = C(x, z)a, the condition is that 
the image of the measure dz on M+ by the function C(x, z) be the measure C(q(3;)) ^i+I(j;) 
and this yields the function 



1 



/ aix)z -a{x)\ °("' 

Our hypotheses on the bottom structure are fulfilled, we may choose 7[j,i*] = A 1)/ 
on Mf^. If we suppose 

a of class n Lip and < Ai ^ a(a;) A2 < 2 (22) 

the assumptions Rla) and Rib) are fulfilled and the condition Rlc) has been proved by 
Hiraba ([H] Remark 3.6 p. 43 et seq.). 

We can conclude that under hypotheses ([22]) there exists a Markov process whose 
transition semi-group admits as generator and this semi-group possesses a density. 



4 Comments. 

The simplest case of SDE with jumps is the case of Levy processes themselves. If y is a 
Levy process with values in the method gives 

m,Y;] = 7[j,r](Ay«) 

aeJT 

where j is the identity. In this case it does not seem that one could do better than the 
Sato criterion [29] (or [30] Thm 27.7). This induces the following natural question: when 
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a Levy measure v is carried by a Lipschitzian curve in M"^ and carries a Dirichlet form 
satisfying <- ^ 7[j, j*](Ays) > a.s. is it necessary such that ^ for some n where 
V = A \)v ? 

When deahng with SDE's driven by Levy processes our approach supposes some reg- 
ularity for the Levy measure of the driver process because of the existence of the bottom 
Dirichlet form and technical condition (HO). It does not need this Levy measure possesses 
a density (cf [7] §2.3 and Example 2 above). Then the existence of density for the solution 
is obtained under weaker hypotheses than those of Leandre (cf [21] [22]) because there is 
no growth condition for the Levy measure near the origin, and also because we do not need 
that the measures on M from which the Levy measure is a sum of images, have C^-densities. 
A similar remark may be done when comparing our hypotheses with those of Ishikawa and 
Kunita [18] who suppose non degeneracy of the Levy measure of the driver Levy process. 
As said in the introduction, the main advantage of the Dirichlet forms method is to allow 
only n Li-p coefficients. Now these authors obtain also smoothness results that we do 
not discuss in the present paper. 

In another work we are studying the extension of our arguments to smoothness results 
thanks to the fact that the gradient defined by formula dl} may be easily iterated. 

Levy processes and random Poisson measures do possess strong regularizing properties 
due to the fact that the jumps are independent of the strict past (cf [7] examples 5.1 to 5.3). 
This phenomenon has been deepened by Fournier and Giet [H] who obtained density for 
the solution of an SDE driven by a Levy process supposing only an absolutely continuous 
Levy measure and without using the Malliavin calculus. Even if their hypotheses on the 
coefficients are slightly stronger than ours, this shows that the Malliavin calculus in the 
spirit of [2] which seems to be the most powerful for this aim, have to be crossed with 
other techniques in order to capture all regularizing properties of Levy processes. 

In this perspective, we believe that the simplification brought by the lent particle 
method gives a tool easier to adapt with various arguments. 
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